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abstract 


An  analytical  model  and  computer  program  has  been  developed  to  analyze 
supersonic,  laminar  separating  flew  fields  over  range.  The  important 
features  of  the  model  are  its  ability  to  include  the  tw  family  nature 
of  the  supersonic  flew  field,  the  matching  of  supersonic  and  subsonic 
profiles  in  the  transonic  region  (which  aids  in  the  selection  of  the  solu- 
tion for  a particular  problem  when  two  exist,  thus  eliminating  the  problem 
of  branching  solutions),  the  inclusion  of  the  normal  momentum  equation 
(and  thus,  normal  pressure  gradients)  throughout  the  flow  field,  and  the 
solution  of  the  separation  problem  by  a marching  technique.  The  ellipti- 
city  of  separated  flow  fields  is  accounted  for  by  casting  the  problem  into 
one  of  the  inverse  nature.  Correlations  with  data  are  utilized  to  determine 
upstream  influence  properties.  Results  for  pressure  distribution  and  heat 
transffer  are  presented  and  compared  with  experimental  data.  The  applica- 
bility/ of  the  model  to  more  general  flow  fields  and  geometries  is  assessed. 
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1 _ INTRODUCTION 

The  inportance  of  the  properties  of  the  boundary  layer  and  its  interaction 
with  the  external  flow  field  far  flews  over  compression  surfaces  is  well 
known  and  the  attendant  role  of  these  effects  on  determining  aerodynamic 
characteristics  of  wing- flap  configurations  and  engine-inlet  ramps  have 
been  previously  assessed.  The  importance  of  the  viscous- invisc id  inter- 
action is  magnified  for  high  Mach  number,  lew  Reynolds  number  flow,  and 
the  design  of  compression  surfaces  must  take  into  account  such  inter- 
action phenomena.  Thus,  the  problen  of  the  interaction  of  supersonic, 
laminar  flow  with  a compression  surface,  and  the  resulting  separated 
flow  field  has  been  the  subject  of  numerous  previous  experimental 
(i.e.  Ref.  1-5)  and  analytical  investigations. 

The  flow  properties  of  shock-boundary  layer  interaction  caused  either 
by  a deflected  flap  or  by  a shock  impinging  on  a boundary  layer  have 
been  shown  to  be  similar  in  nature  and  correlative.  The  concept  of  free 
interaction6  allows  for  the  atployment  of  empirical  relations  from  a single 
set  of  free  interaction  experimental  data  to  predict  the  occur enoe  and  the 
extent  of  a separated  zone,  as  well  as  the  pressure  distribution  in  the 
region  of  separation.  It  has  been  shown  experimentally7  that  in  the  region 
of  free  interaction  the  heat  transfer  profiles  are  independent  of  the  mechanism 
by  which  they  were  induced. 

In  order  to  totally  quantify  the  entire  flow  field,  theoretical  calculations 
have  had  to  deal  with  the  basic  ellipticity  of  the  problem.  The  positive 
pressure  disturbance  caused  by  the  shock  (or  by  the  downstream  geo- 
metry) propagates  upstream  through  the  subsonic  portion  of  the  boundary 
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layer  causing  the  boundary  layer  to  separate.  Because  of  the  basic  elliptic 
nature,  it  has  been  inpcrsible  until  recently  to  arplay  a forward  marching 
scheme  without  same  sort  of  iteration  an  the  initial  conditions. 

Marry  analytical  efforts  over  the  past  fifteen  years  have  been  based  on 
integral  techniques  8-H  have  used  a moment  of  conservation  equations-^  »i3 
to  yield  an  undetermined  parameter  which  is  adjusted  to  conform  with  the  up- 
stream influence  effect.  The  method  of  Ref.  13  which  uses  a one  parameter 
family  of  profiles,  gives  good  agreement  with  experimental  pressure  data  for 
adiabatic  and  moderately  cooled  walls,  however,  for  cold  walls  the  method 
has  proved  inadequate.  Also,  the  heat  transfer  results  obtained  by  this 
method  for  all  wall  temperature  ratios  have  shewn  poor  quantitative  agree- 
ment with  experiments.  Holden^  has  extended  the  method  of  Ref.  13  to 
allcw  for  separate  parameters  for  the  velocity  and  enthalpy  profiles 
improving  agreement  with  experiment  for  the  highly  cooled  cases.  All  of 
the  integral  techniques  have  neglected  the  possible  effects  of  a normal 
pressure  gradient. 

More  recent  studies^ #15,16  have  attempted  to  use  finite  difference  tech- 
niques coupled  to  an  interaction  pressure  law  for  the  inviscid  flow. 

These  solutions  have  shown  a tendency  to  have  an  infinite  number  of  branch- 
ing solutions  possible  17,18^  These  branched  solutions  were  explained  as 
being  caused  by  various  types  of  downstream  geometries  causing  different 
rates  and  magnitues  of  upstream  influence.  DwsyeilG  investigated  a variety 
of  branch  selections  aixl  determined  that  to  control  downstream  reattachment 
and  relaxation  to  downstream  conditions,  it  was  most  convenient  to  preselect 
the  second  derivative  of  displacement  thickness  at  the  initial  station. 


This  parameter  could  then  be  adjusted  by  an  iterative  procedure  until 
the  proper  downstream  conditions  were  recovered.  All  such  methods  have 
also  neglected  the  possible  effects  of  normal  pressure  gradients  which 
may  be  large  in  the  region  of  separation  and  the  effects  of  a two- family 
supersonic  inviscid  flow. 

A recent  investigation10  has  attainted  to  develop  an  algorithm  which 
eliminates  the  necessity  of  iteration  by  using  the  downstream  boundary 
conditions  and  recasting  the  interaction  equations  into  a tune  dependent 
set.  The  steady  state  conditions  then  corresponds  to  the  physical  prob- 
lan  of  interest.  All  previous  studies,  whether  utilizing  finite  differ- 
ence or  integral  techniques  have  had  to  assume  a pressure-streamline 
deflection  relationship  at  the  boundary  layer  edge  (either  a linear  depen- 
dence or  a tangent  wedge  approximation) , to  determine  interaction  effects 
and  thus  have  of  necessity,  neglected  the  effects  of  embedded  shocks, 
and  the  attendant  effects  of  downrunning  waves  which  can  alter  boundary 
layer  properties.  In  addition,  whereas  a linear  pressure-deflection 
relationship  applies  for  weak  interaction  and  the  tangent  wedge  relation- 
ship approximates  strong  interaction,  it  is  always  a problem  to  decide 
for  a particular  problem  which  relationship  to  use, especial  ly  in  the 
overlap  domain  (i.e.  X~0  (1)). 

Previous  work  at  New  York  University20'21'22  has  attempted  to  establish  a 
unified  approach  for  the  solution  of  expanding  and  oppressing  flews 
where  the  effects  of  interaction  are  important.  The  investigations  have 
sought  to  demonstrate  the  validity  of  the  model  developed  (which  will  be 
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described  in  the  next  section)  and  the  general  applicability  of  the  numer- 
ical analysis  which  has  been  produced.  The  important  features  of  the 
model  are:  its  ability  to  include  the  two  family  feature  of  the  super- 
sonic field  (which  all  other  previous  analyses  have  neglected) , the  match- 
ing of  supersonic  and  subsonic  profiles  in  the  transonic  region  (which  aids 
in  the  selection  of  the  solution  for  a particular  problem  when  two  exist, 
thus  eliminating  the  problem  of  branching  solutions)  and  the  inclusion  of 
the  normal  momentum  equation  (and  thus  nomal  pressure  gradients)  in  the  super- 
sonic field  (by  the  use  of  the  viscous  characteristic  equations)  and  in  the 
highly  viscous  subsonic  zone  (solved  by  an  irrplicit  finite  difference  scheme) . 
It  has  thus  been  asserted^  that  interaction  problems  can  be  well  posed  as 
initial  value  problems  within  the  context  of  this  model. 

The  ellipticity  of  separating  and  reattaching  flews  over  ramps  adds  on  addi- 
tional parameter  to  the  general  description  of  the  model,  which  is  based 
on  a marching  scheme.  The  effects  of  strong  upstream  influence  must  be 
accounted  for  and  the  following  sections  delineate  hew  this  is  handled. 
Following  a section  describing  the  model  and  its  special  details  necessary 
for  separated  flow  fields,  a section  on  the  equations  of  motion  in  each  of 
the  regions  of  the  problem  is  presented.  Then,  results  of  the  analysis  for 
two  separated  ramp  flow  fields  are  described.  The  results  are  compared  with 
experiments,  and  conclusions  of  the  model's  applicability  to  separated 
flow  fields  is  presented. 


' 


I 
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II.  DESCRIPTION  OF  THE  PROBLEM  AND  THE  ANALYTICAL  MODEL 

Ihe  basic  problan  which  is  studied  is  a supersonic  or  hypersonic  flow  over 
a wing-flap  configuration.  As  an  example  of  this  phencnvanon  consider  a 
flew  over  a flat  plate-wedge  configuration  shown  schematically  in  Fig.  1. 

Ihe  start  of  the  free  interaction  region  occurs  at  the  station  (Xo)  where 
the  wall  properties  begin  to  depart  from  flat  plate  values.  In  this 
region,  an  adverse  pressure  gradient  is  aonmunicated  through  the  subsonic 
portion  of  the  boundary  layer,  causing  the  wall  pressure  coefficient  to 
rise  while  simultaneously,  the  skin  friction  and  heat  transfer  decrease. 

Outgoing  compression  waves  ber  go  feed  upward  into  the  freestream. 

For  supersonic  Mach  numbers,  this  compression  fan  coalesces  into  a shock 
sate  distance  away  from  the  wall.  In  the  hypersonic  case,  the  shock 
forms  almost  at  the  boundary  layer  edge  (or  near  the  subsonic- supersonic 
viscous  flew  interface) , its  slope  approximating  the  flow  deflection  of 
the  boundary  layer  edge  itself.  In  both  cases  the  flow  downstream  of 
the  coalescence  may  be  dominated  by  a two  family  characteristic  flow  field 
as  dewnrunning  waves  from  the  shock  may  significantly  alter  the  pressure 
distribution  given  by  a Prandtl-Meyer  relation.  Hie  separation  point  (Xs) 
is  defined  by  the  point  of  vanishing  skin  friction.  A dividing  stream- 
line can  be  identified  which  forms  the  boundary  between  positive  and 
negative  values  of  the  stream  function.  A second  compression  fan  and 
shock  are  found  due  to  the  reattachment  compression. 

Figure  2 depicts  schanatically  the  wall  pressure,  heat  transfer  and  skin 
friction  distributions.  In  the  region  of  free  interaction  the  pressure 
gradient  is  approximately  constant  and  then  decreases  gradually  until  a 

pressure  plateau  is  reached.  In  the  region  of  the  flap  the  pressure  again  i 
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begins  to  rise,  reaching  a maximum  and  then  decreasing  asymptotically  to 
the  two-dimensional  wedge  value.  The  skin  friction  after  the  hinge 
increases  fran  a negative  value,  passes  through  zero  at  reattachment 
and  increases  again  to  the  flat  plate  value.  The  heat  transfer  distri- 
bution reaches  its  minimum  in  the  region  of  the  hinge  and  then  increases 
as  shewn  in  the  figure. 

The  analytical  model  and  method  of  analysis  developed  in^O, 21,22  for  ejq?and- 
ing  flow  fields  has  been  modified  to  allow  for  the  analysis  of  the  above 
problem.  The  method  of  analyses  is  as  follows:  In  the  region  prior  to 

the  effect  of  upstream  influence  the  method  described  in  Ref.  21  is  used 
unaltered.  The  supersonic  portion  of  the  flow  field  is  solved  by  charac- 
teristics which  have  been  modified  to  include  viscous  transports  normal  to 
the  streamlines.  The  inclusion  of  viscosity  in  the  characteristic  rela- 
tions allows  for  an  inner-outer  matching  to  occur  in  the  transonic  region. 

The  equations  for  the  inner  region  (the  highly  viscous  subsonic  flow  near 
the  wall)  are  represented  by  the  boundary  layer  equations  in  coordinates 
measured  along  and  normal  to  the  surface.  The  dictate  from  the  outer 
flow  field  that  normal  pressure  gradients  may  occur,  coupled  with  the  effect 
of  the  thickening  boundary  layer  (enhanced  for  moderate  Reynolds  numbers 
and/or  the  onset  of  separation)  leads  to  the  inclusion  of  normal  pressure 
gradient  effects  in  the  inner  flow  field.  Thus,  the  second  momentum  equa- 
tion is  retained  in  the  inner  region.  The  inner  equations  are  solved  by  an 
implicit  finite  difference  Crank -Nicholson  scheme.  The  matching  between  inner 
and  outer  solutions  is  acccnplished  by  guessing  a streamline  deflection  at 
the  matching  point  of  the  next  step  (whose  stepsize  is  dictated  by  the  charac- 
teristic intersection  distance)  and  determining  the  pressure  at  the  tran- 
sonic intersection  from  the  downrunning  characteristic  relation.  The  inner 


region  is  then  ccmputed  and  the  condition  that  at  the  wall  the  stream 
function  is  zero  is  tested.  An  iteration  procedure  based  on  changing  the 
streamline  deflection  is  continued  until  the  surface  condition  (l^all  = 0) 
is  met. 

The  interaction  program  is  then  applied  to  the  separation  problem  in  an 
inverse  manner.  At  a selected  station  along  the  flat  plate,  a selected 
adverse  pressure  gradient  is  applied  at  the  wall  (The  method  for  the  selec- 
tion of  Xo  and  is  explained  in  Section  IV. ) Consistent  with  the  thus 
known  wall  pressure,  the  iteration  procedure  is  dropped.  This  is  equiva- 
lent to  allowing  for  mass  injection  at  the  surface  (since  ’i'wall^O) . It 
has  previously  been  shewn ^ that  there  is  a close  correspondence  between 
the  upstream  influence  of  a downstream  flap  and  a boundary  layer  with 
injection.  This  procedure  is  continued  until  flow  separation.  At  the 
separation  point,  the  complete  iteration  procedure  is  again  reverted  to 
(it  is  no  longer  necessary  to  specify  the  adverse  pressure  gradient)  and 
the  syston  is  solved  consistent  with  ^wall  = 0.  An  inportant  result  is 
that  the  pressure  gradient  which  is  determined  in  the  course  of  the  march- 
ing is  entirely  consistent  with  the  imposed  adverse  pressure  gradient 
prior  to  separation.  It  is  the  assertion  here  that  the  pressure  gradient 
is  continuous  after  separation  due  to  the  fact  that  the  supersonic  charac- 
teristic field  in  the  vicinity  of  the  matching  has  acoortmodated  to  the 
effect  of  streamline  curvature  and  thus  the  dewnrunning  characteristic 
(which  determines  the  pressure)  is  consistent  with  the  previous  pressure 
distribution.  An  embedded  shock  routine  has  been  added  to  the  main  pro- 
gram to  allow  for  characteristic  coalescence  when  it  occurs  (See  Fig.  1). 
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It  has  been  shown  14,16  that  instabilities  can  arise  in  the  numerical 
solution  of  the  boundary  layer  equations  by  a marching  technique  in 
regions  of  reverse  flew.  This  arises  frem  the  fact  that  the  convective 
terms  puux  and  PuTx  in  the  momentum  and  energy  equations  become  nega- 
tive. These  instabilities  also  arose  in  the  present  analysis  in  the 
region  of  separation.  It>  overcome  this  effect,  the  same  approach  as 
has  been  used  in  previous  investigations  was  employed,  i.e,  the  utili- 
zation of  the  absolute  value  of  u in  the  convective  terms 
(p  J u | ux  and  p |u|Tx) . 

At  a certain  previously  unknown  station,  the  program  is  no  longer  able 
to  converge  on  a solution  (i.e.,  pass  the  mass  flow).  At  this  station 
the  wall  is  turned.  The  position  of  and  the  angle  of  the  turn  which 
were  previously  unknown  are  the  two  parameters  which  are  now  determined, 
thus  rendering  the  program  one  of  an  inverse  nature.  Hius,  the  ellipticity 
of  the  flow  field  is  overcome  by  determining  the  surface  geometry  which 
caused  the  particular  upstream  separating  pressure  gradient  (and  its  posi- 
tion) as  a consequence  of  the  solution.  The  wall  turning  is  accomplished 
within  the  program  by  turning  the  wall  through  a snail  A0  in  each  step  Ax 
until  it  is  no  longer  necessary  to  change  to  pass  the  mass  flow.  Thus, 
the  angle  of  wall  turning  is  determined  and  the  program  is  continued 
through  reattachment. 

Thus,  the  flew  field  has  been  separated  into  six  distinct  regions  of  ooirpu- 
tation  (Fig.  3)  . 

I.  Region  prior  to  upstream  influence 
II.  Region  of  prescribed  adverse  pressure  gradient 


Region  of  separated  flow  where  convective  terms  are  utilized  with 
absolute  value  of  u in  the  negative  region 
Region  of  wall  turning  (dictated  by  the  inability  to  have  a 
converged  solution) 

Region  approaching  reattachment  along  wedge  (solved  in  same  way 
as  region  III) 

Region  of  attached  flow  relaxing  to  Blasius  solution 


III.  EQUATIONS  OF  MOTION 


The  equations  of  notion  used  in  the  distinct  regions  of  the  flow  field 
have  been  presented  previously  2 0-2 2 ^ in  the  supersonic  portion  of  the 
flow  field  viscous  transports  normal  to  the  streamlines  cure  included  in 
the  characteristic  yielding  the  following  two-dimensional  equations  for 
the  determination  of  the  properties: 


curvature  terms  neglected) 


(P  u)x  + (pv)f  = 0 


p ,'uu  x + VU^  ) +PX  ~ 


P/  UV  x 


+ 


vv- 

J — iL 

/ Rea, 


+ Py  i^Ffeoo  — 0 


(3) 

(4) 

(5) 


p (u  Tx  + v Ty  ) - Ea,  (u  px  + ^Py)  = S2 ' (6) 

where  = (V  ^y 

S'  = — (pT-)-  + EooP  (Uy)  2 

2 Pr  y'y  y 

The  solution  of  the  systen  of  equation  and  method  of  iteration  valid  in 
each  region  of  the  flow  field  has  been  described  in  the  previous  section. 
One  additional  statement  must  be  made  with  respect  to  the  solution  of  the 
continuity  equation  for  v in  the  subsonic  zone.  The  matching  of  solutions 
in  the  transonic  zone  causes  one  to  be  ccrpelled  to  use  a small  stepsize 
Csince  in  the  transonic  zone  characteristics  are  close  to  vertical) . This 
snail  stepsize  has  caused  an  instability  in  the  determination  of  v from 

, < dip 

continuity  when  it  is  computed  from  the  equation  pv  = - -jjj  . This  has 
been  overcome  consistently  for  both  expanding  and  compressing  flow  fields 
by  utilizing  the  modified  equation  Pv  = ~ 355  • where  AX  xs  taken  to  be  an 

the  order  of  ten  AX  and  then  decreased  during  the  course  of  the  marching. 

In  the  following  section  results  of  the  computer  program  are  presented  and 
oonparisons  with  experiments  are  shown  to  demonstrate  the  validity  of  the 
model. 
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IV  . RESULTS  AND  CONCLUSIONS 

The  basic  ellipticity  of  the  problan  necessitates  a determination  of  two 

parameters,  the  point  of  demarkation  frcm  the  flat  plate  results  (Xo)  and 

the  pressure  gradient  frcm  there  to  the  separation  point  f , in  order 

xaxJ  Xo 

to  proceed  with  the  calculation.  The  problem  is  thus  cast  as  one  of  an 
inverse  nature,  i.e. , select  and  Xo  and  and  find  the  plate  length  L and 
wall  turn  6^  which  brought  about  the  upstream  effect.  In  order  to  calcu- 
late the  flow  for  a particular  problan  it  is  thus  necessary  to  iterate  on  the 
initial  conditions  until  the  tvo  parameters  are  consistent  with  the  flow  of 
interest. 

In  order  to  somewhat  overcome  this  complexity,  and  large  increase  in  computer 
time  per  case,  an  analysis  of  available  data  was  undertaken  to  correlate 
results  for  these  parameters  over  a wide  range  of  the  interaction  parameter 
X. 

24  25  26 

Based  on  the  suggestions  of  ' , Cheng  D found  that  the  inportant  parameter 

2. 

governing  the  similitude  of  the  flow  at  separation  is  (M»0  incip)/x£  where 

Gincip  is  the  ramp  angle  at  incipient  separation  and  xe  = (.664  + 1.73'IW/TC^X> 

• • • 

This  effect  was  correlated  over  a wide  range  from  weak  to  strong  interaction.  It 
was  also  shown  that  the  upstream  influence  in  an  attached  flow  was  inversely  pro- 
portional to  Mo/  /Rbl- 


Based  on  these  conclusions,  the  data  of  separated  flew  fields  was  investigated 
to  determine  the  oorre lativity  of  and ^ as  a function  of  Xc*  It 

^ ly*  Q X VA/t 

was  found  that  for  0^  > ©incip  (separated  flows),  the  parameter  — 

governs  the  similitude  of  the  flow  and  that  in  the  region  e 

1/2 

XO  < X < Xe  also  varies  as  x.  , These  results  are  presented  in  Figures 
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4 and  5,  which  present  data  over  a wide  range  of  the  interaction  param- 
eter. 

These  correlations  thus  allow  for  the  selection  of  the  two  parameters  Xo  and 

dp  for  any  flow  condition  and  ramp  angle  with  reasonable  accuracy.  Hie  iter- 
dx 

at ion  procedure  to  determine  the  exact  values  for  a particular  geometry  is 
thus  reduced  considerably  and  the  inverse  characteristic  of  the  model  no 
longer  produces  a major  drawback. 

The  resilts  of  two  calculations  corresponding  to  the  data  of  Ref.  2 are  pre- 
sented here.  For  these  cases,  the  correlations  utilized  allowed  that  no 
iteration  on  initial  conditions  were  needed  to  find  the  flow  fields  consis- 
tent with  the  cases  to  be  reproduced.  The  results  for  wall  pressure  and 
heat  transfer,  presented  in  Figures  6-8,  indicate  that  the  model  reproduces 
the  flow  phenomena  with  good  accuracy. 

It  is  thus  the  contention  that  separated  flow  fields  in  the  supersonic  and 
hypersonic  regimes  can  be  computed  utilizing  the  interaction  model  presented 
herein,  which  computes  the  entire  flow  at  each  step  of  a inarching  technique 
without  the  necessity  of  specifying  a pressure-deflection  relationship  and 
allowing  for  anbedded  shocks  and  the  two  family  nature  of  the  flow  field. 

The  inverse  nature  of  the  analysis  is  somewhat  overcome  for  ramps  by  utilizing 
the  correlations  developed.  This,  of  course,  would  not  be  the  case  for  more 
complex  geometries.  The  major  drawback  of  this  model  is  the  necessity  for 
a given  flow  condition,  to  have  an  available  correlation  for  the  determination 
of  the  upstream  influence  effect.  Therefore,  while  in  principle  the  model 
has  general  applicability  for  the  solution  of  laminar  separated  flowfields 
in  the  supersonic  and  hypersonic  regime,  in  practice,  complex  geometric  con- 
figurations (including  those  which  yield  three-dimensional  flowfields)  may  be 
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mare  easily  accomodated  within  the  scope  of  fully  elliptic  programs. 

With  respect  to  turbulent  flows,  one  could  foresee  a similar  situation. 
Ample  experimental  data  is  available  for  turbulent  flows  over  ranps  to 
develop  a correlation  similar  to  the  ones  presented  in  Figures  4 and  5. 

Che  would  expect  that  the  correlations  would  display  a different  power 
law  than  in  the  laminar  case,  perhaps  yielding  solutions  which  vary  as  X-8. 
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SCHEMATIC  REPRESENTATION  OF  SUPERSONIC 


FIGURE  3.  DISTINCT  REGIONS  OF  COMPUEATICN 


FIGURE  4.  CORRELATION  OF  DATA  IN  TERMS  OF  INTERACTION  PARAMETER  TO  DETERMINE  SEPARATION  Xq/L 


Re,  xIO 


FIGURE  7.  PRESSURE  DISTRIBUTION  IN  A SEPARATED  FLOW  FT  LID  AND  COMPARISON  WITH  DATA  OF  REF. 
= 7.5,  Fte  = 2.9  x 105/ft.,  TW/To  = .19,  0„  = 20°) 


FIGURE  8.  HEAT  TRANSFER  DISTRIBUTION  IN  A SEPARATED  FLOW  FIELD  AND  COMPARISON  WITH  DATA  OF  REF. 


REFERENCES 


1.  Miller,  D . Hynan,  R.  and  Childs,  M , "Mach  8 to  22  Studies  of  Flow 
Separations  IXue  to  Deflected  Control  Surfaces",  AIAA  Journal,  Vol  2, 

No.  2,  Feb.  1964 

2.  Johnson,  C B , "A  Theoretical  and  Experimental  Study  at  Mach  8 of  Flow 
Separation  of  a Flat  Plate  With  Deflected  Trailing  Edge",  M S -Thesis, 

V PI.,  Dept,  of  Aerospace  Engineering,  Aug.  1966 

3.  Lewis,  J .Kubota,  T and  Lees,  L , "Experimental  Investigation  of  Super- 
sonic, laminar  Two-Dimensional  Boundary  Layer  Separation  in  a Compression 
Comer  With  and  Without  Cooling",  AIAA  Journal,  Vol.  6,  No.  1,  Jan.  1968 

4.  Sfeir,  A , "Supersonic  Laminar  Boundary  Layer  Separation  Near  a Compression 
Comer",  Dept.  No.  AS 6-6,  U.  Cal.  at  Berkeley,  March  1969 

5.  Holden,  M ."Leading  Edge  Bluntness  and  Boundary  Layer  Displacement  Effects 
on  Attached  and  Separated  Laminar  Boundary  Layers  in  a Compression  Comer" , 
CAL  -Report,  1969 

6.  Erdos,  J .and  Pal  lone,  A . "Shock- Boundary  Layer  Interaction  and  Flow 
Separation" , Proc.  1962  Heat  transfer  and  Fluid  Mechanics  Institute, 

Stamford  U.  Press,  1962 

7.  Gadd,  G E , "Boundary  layer  Separation  in  the  Presence  of  Heat  Transfer", 
AGARD  Report  No.  280,  April  1960 

8.  Curie,  N.  "The  Effects  of  Heat  Transfer  on  Iaminar-Eoundary- Layer  Separation 
in  Supersonic  Flow",  The  Aeronautical  Quarterly,  Vol.  XII,  pp.  309-336, 
November,  1961 

9.  Savage,  S.B. , "The  Effect  of  Heat  Transfer  on  Separation  of  laminar  Ccrpres- 
sible  Boundary  layers",  CALCIT  Separated  Flows  Project,  Technical  Report 
No.  2,  June  1,  1962 

10.  Glick,  Herbert  S. , "Modified  Crocco-Lees  Mixing  Theory  for  Supersonic 
Separated  and  Reattaching  Flows",  GAICIT  Hypersonic  Research  Project 
Memorandum  No.  53,  May  2,  1960 

11.  Abbott,  D.  E , Holt,  M.,  and  Nielsen,  J.N  , "Investigation  of  Hypersonic 
Flow  Separation  and  its  Effect  on  Aerodynamic  Control  Characteristics", 
ASD-TDR-62-963,  Nov.  1962 

12.  Lees,  L. , and  Reeves,  L. , "Supersonic  Separated  and  Reattaching  Laminar  Flows: 
I .General  Theory  Application  to  Adiabatic  Boundary  Layer-Shock  Wave  Inter- 
actions", AIAA  Journal,  Vol.  2,  No.  11,  Nov.  1964,  pp.  1907-1920 

13.  Holden,  M.S. , "An  Analytical  Study  of  Separated  Flows  Induced  by  Shock  Wave 
Boundary  Layer  Interaction",  CAL  Report  No.  AI-1972-A-3,  Dec.  1965 

14.  Rehner,  T.A  .and  Flugge-Lotz,  I , "The  Interaction  of  a Shock  Wave  with  a 
Laminar  Boundary  Layer",  Int.  Journal  of  Non-Linear  Mechanics,  Vol.  3,  1966 


% 


-23- 


A 


15.  Werle,  M.  J.,  FOlak,  A.  and  Bertke,  S.D. , "Supersonic  Boundary  Layer 
Separation  and  Reattaclment  Finite  Difference  Solutions" , Kept.  AFL 
72-12-1,  Dept,  of  Aeronautical  Ihgineering,  University  of  Cincinnati, 
January  1973 

16.  Dwayer,  D.L. , "Supersonic  and  Hypersonic  Two-Dimensional  Laminar  Flow 
Over  a expression  Comer" , Proc.  AIAA  Ccrputationa  1 Fluid  Dynamics 
Conference,  July  1973 

17.  Baum,  E , "An  Interaction  Model  of  a Supersonic  Laminar  Boundary 
Layer  on  Sharp  and  Rounded  Backward  Facing  Steps",  AIAA  Journal, 

Vol.  6,  No.  3,  March  1968 

18.  Werle,  M.J  -and  Bertke,  S.D. , "A  Finite  Difference  Method  for  Boundary 
Layers  with  Reverse  Flew",  AIAA  Journal,  Vol.  lo.  No.  9 

19.  Werle,  M.J.  and  Vatsa,  V.N. , "Numerical  Solution  of  Interacting  Super- 
sonic Boundary  layer  Flows  Including  Separation  Effects",  A. R.L. -73-0162, 
Decenber  1973 

20.  Ferri,  A.  and  Dash  S.,  "Viscous  Flew  at  High  Mach  Numbers  With  Pressure 
Gradients",  published  in  "Viscous  Interaction  Phenomena  in  Supersonic 
and  Hypersonic  Flow",  University  of  Dayton  Press,  1969 

21.  Miller,  G. , "The  Mathenatic  Formulation  of  Viscous- Invisc id  Problems  in 
Supersonic  Flow",  AIAA  Journal,  Vol.  11,  No.  7,  July  1973 

22.  Miller,  G.  and  Srowkcwski,  A.,  "Model  for  the  Solution  of  Supersonic  Vis- 
cous Interaction  Effects  on  Blunted  Cones",  to  appear  in  AIAA  Journal, 
Feb.  1975 

23.  Behrens,  W. , "Separation  of  a Supersonic  Turbulent  Boundary  Layer  By  a 
Forward  Facing  Step",  AIAA  Paper  71-127,  January,  1971 

24.  Holden,  M.S. , "Theoretical  and  Experimental  Studies  of  Separated  Flows 
Induced  by  Shock  Wave- Boundary  Layer  Interaction",  A3ARD  Specialists 
Meeting  on  Separated  Flews,  Conference  Proceedings  No.  4,  Part  I, 

May  1966 

25.  Needham,  D.A. , "laminar  Separation  in  Hypersonic  Flow",  University 
of  London  Ph.D.  Thesis  (August  1965) 

26.  Cheng,  H.K. , Hall,  J.G. , Golian,  T.C. , and  Hertzberg,  A.,  "Boundary 
Layer  Displacement  and  Leading-Edge  Bluntness  Effects  in  High-Tenper- 
ature  Hypersonic  Flew",  J.Aero/Space  Sci.  Vol.  28,  No.  5,  May  1961 


-24- 


1 


